Applications of complementarity plot in error detection and structure validation of proteins.

The complementarity plot (CP) is based on packing and electrostatics of amino acid residues buried within globular proteins and is a sensitive indicator of the harmony or disharmony of interior residues with regard to short and long range forces sustaining the native fold. As a structure validation tool, it has already been reported to be effective in detecting erroneous side-chain torsions in obsoleted structures. The current study describes the design of several local and global scores based on CP and surveys their utilities in discriminating between obsolete structures and their corresponding upgraded counterparts, detection of wrong rotamer assignments and in identifying packing anomalies. CPs are especially effective in the detection of low-intensity errors (in main-chain geometrical parameters) diffused over the entire polypeptide chain. The methodology is also used to confirm the integral role played by strategic deviations (in main-chain geometrical parameters) in maintaining fold integrity, as reversal to their corresponding ideal values (either unimodal or conformation dependent) lead to large-scale structural distortions. A special feature of this validation tool is to signal unbalanced partial charges within protein interiors. The application of CP in protein homology modeling and protein design is also demonstrated.

Effect of temperature on DNA double helix: An insight from molecular dynamics simulation.

The three-dimensional structure of DNA contains various sequence-dependent structural information, which control many cellular processes in life, such as replication, transcription, DNA repair, etc. For the above functions, DNA double helices need to unwind or melt locally, which is different from terminal melting, as often seen in molecular dynamics (MD) simulations or even in many DNA crystal structures. We have carried out detailed MD simulations of DNA double helices of regular oligonucleotide fragments as well as in polymeric constructs with water and charge-neutralizing counter-ions at several different temperatures. We wanted to eliminate the end-effect or terminal melting propensity by employing MD simulation of DNA oligonucleotides in such a manner that gives rise to properties of polymeric DNA of infinite length. The polymeric construct is expected to allow us to see local melting at elevated temperatures. Comparative structural analysis of oligonucleotides and its corresponding virtual polymer at various temperatures ranging from 300 K to 400 K is discussed. The general behaviour, such as volume expansion coefficients of both the simulations show high similarity, indicating polymeric construct, does not give many artificial constraints. Local melting of a polymer, even at elevated temperature, may need a high nucleation energy that was not available in the short (7 ns) simulations. We expected to observe such nucleation followed by cooperative melting of the polymers in longer MD runs. Such simulations of different polymeric sequences would facilitate us to predict probable melting origins in a polymeric DNA.

Theoretical analysis of noncanonical base pairing interactions in RNA molecules.

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.

Energy barriers and rates of tautomeric transitions in DNA bases: ab initio quantum chemical study.

Tautomeric transitions of DNA bases are proton transfer reactions, which are important in biology. These reactions are involved in spontaneous point mutations of the genetic material. In the present study, intrinsic reaction coordinates (IRC) analyses through ab initio quantum chemical calculations have been carried out for the individual DNA bases A, T, G, C and also A:T and G:C base pairs to estimate the kinetic and thermodynamic barriers using MP2/6-31G** method for tautomeric transitions. Relatively higher values of kinetic barriers (about 50-60 kcal/mol) have been observed for the single bases, indicating that tautomeric alterations of isolated single bases are quite unlikely. On the other hand, relatively lower values of the kinetic barriers (about 20-25 kcal/mol) for the DNA base pairs A:T and G:C clearly suggest that the tautomeric shifts are much more favorable in DNA base pairs than in isolated single bases. The unusual base pairing A':C, T':G, C':A or G':T in the daughter DNA molecule, resulting from a parent DNA molecule with tautomeric shifts, is found to be stable enough to result in a mutation. The transition rate constants for the single DNA bases in addition to the base pairs are also calculated by computing the free energy differences between the transition states and the reactants.